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Lattice simulations employing reweighting and Taylor expansion techniques have predicted a 
(ji , T) -phase diagram according to general expectations, with an analytic quark-hadron crossover 
at jl = turning into a first order transition at some critical chemical potential By contrast, re- 
cent simulations using imgainary ji followed by analytic continuation obtained a critical structure 
in the {m u d,m s ,T,jj,} parameter space favouring the absence of a critical point and first order line. 
I review the evidence for the latter scenario, arguing that the various raw data are not inconsistent 
with each other. Rather, the discrepancy appears when attempting to extract continuum results 
from the coarse (Nt = 4) lattices simulated so far, and can be explained by cut-off effects. New (as 
yet unpublished) data are presented, which for Nf — 3 and on N t = 4 confirm the scenario without 
a critical point. Moreover, simulations on finer Nt = 6 lattices show that even if there is a critical 
point, continuum extrapolation moves it to significantly larger values of jj.e than anticipated on 
coarse lattices. 
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Figure 1: Left: Expected phase diagram for physical QCD. Right: The prediction from a recent lattice in- 
vestigation [^3J. The entire phase boundary between plasma and confined phase corresponds to a crossover. 

1. Introduction 

One of the most fascinating and challenging tasks in the study of QCD under extreme condi- 
tions remains to elucidate its phase diagram as a function of temperature and baryon density. QCD 
is a strongly coupled theory, and it has recently become clear that certain strong coupling features 
persist also in the quark gluon plasma regime up to the temperatures that are currently accessible 
experimentally. Hence, perturbation theory fails us for this problem, and simulations of lattice 
QCD are the only known tool by which we may hope to ultimately come to reliable theoretical 
predictions for QCD. Unfortunately, until 2001 the so-called "sign-problem" of lattice QCD has 
prohibited simulations at finite baryon density. 

Despite this lack of reliable calculations, it is generally believed that the phase diagram quali- 
tatively looks as in Fig. [I] (left). Lattice simulations have been operating along the temperature axis 
at jJ. = for a long time, and results have been consistent with an analytic crossover in this regime. 
However, it was only very recently that these computations could be performed with realistic quark 
masses and on reasonably fine lattices [Jin, so that we now may have confidence in this prediction. 
For T = and finite density, on the other hand, a number of models and perturbative arguments 
[Q, ||] predict a first order phase transition to a superconducting regime at some fi c > 1 GeV. In the 
simplest scenario the phase boundary encountered in this finite density transition and that of the 
finite temperature transition are continuously connected. In this case the first order transition line 
must terminate in a second order end-point, thus leading to the familiar phase diagram, Fig. [j] (left). 
For more detailed arguments and references, see e.g. [§, 

Significant progress in lattice QCD was made since 2001, with several different methods now 
available that circumvent the sign problem, rather than solving it: i) Multi-parameter reweighting 
[§], ii) Taylor expansion around /x = [§] and iii) simulations at imaginary chemical potential, 
either followed by analytic continuation [f7|] or Fourier transformed to the canonical ensemble [j|]. 
It is important to realize that all of these approaches introduce some degree of approximation. 
However, their systematic errors are rather different, thus allowing for powerful crosschecks. All 
methods are found to be reliable as long as /T < 1 , or < 550 MeV, which includes the region of 
interest for heavy ion collisions. Reviews specialized on technical aspects can be found in [||, |T0| ]. 

In particular, a calculation using reweighting methods for Nf = 2 + 1 with physical quark 
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Figure 2: Left: Schematic phase transition behaviour of Nt = 2 + 1 flavor QCD for different choices of 
quark masses (m u j,m s ), at fi = Right: Numerical results for the chiral critical line. Also shown are the 
physical point according to Jll[], and a fit corresponding to a possible tricritical point m'" c ~ 2.8T. 



masses on a coarse lattice (N t = 4, a ~ 0.3 fm) appeared to confirm the expected phase diagram and 
gave a prediction for the critical end point at /if ~ 360(40) MeV [p]. Similarly, an investigation 
of the convergence radius for the Taylor series of the pressure in the Nf = 2 theory with bare quark 
mass m/T = 0.1 yielded a sig nal for the endpoint at H§/T E = 1.1 ±0.2, T E /T c (pL = 0) = 0.95 
[12]. By contrast, a recent study of Nf ■ = 2 + 1,3 using imaginary chemical potentials reached a 
surprising conclusion contradicting the standard scenario: for physical values of the quark masses 
there would be no critical point or first order line at all, the crossover region extends all the way to a 
possible line delimiting the superconducting phase, as in Fig. |] (right) [13]. In the following I will 
go through the evidence for this unorthodox scenario and argue that the raw data of the apparently 
contradicting investigations are actually consistent with each other. The apparent discrepancy arises 
when attempting to extract continuum physics and can be explained by standard cut-off effects. 



2. The nature of the QCD phase transition: qualitative picture 

Rather than fixing QCD to physical quark masses and just varying T,/j., it is instructive to 
consider varying quark masses and the full parameter space {m U:C i,m s ,T ', ju } . The qualitative picture 
for the order of the phase transition at zero density is shown in Fig. || (left). In the limits of zero 
and infinite quark masses (lower left and upper right corners), order parameters corresponding to 
the breaking of a global symmetry can be defined, and one finds numerically at small and large 
quark masses that a first-order transition takes place at a finite temperature T c . On the other hand, 
one observes an analytic crossover at intermediate quark masses. Hence, each corner must be 
surrounded by a region of first-order transition, bounded by a second-order line. 

Now consider the effect of a baryonic chemical potential, /!# = 3/1. As a function of quark 
chemical potential jU, represented vertically in Fig. [| the critical line determined at /I = now 
spans a surface. The standard expectation for the QCD phase diagram is depicted in Fig. || (left). 
The first order region expands as /I is turned on, so that the physical point, initially in the crossover 
region, eventually belongs to the critical surface. At that chemical potential /i£, the transition is 
second order: that is the QCD critical point. Increasing jj. further makes the transition first order. 
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Figure 3: The chiral critical surface in the case of positive (left) and negative (right) curvature. If the 
physical point is in the crossover region for fi = 0, a finite fi phase transition will only arise in the scenario 
(left) with positive curvature, where the first-order region expands with ji. Note that for heavy quarks, the 
first-order region shrinks with fi [|l4|. 



Drawn in the (r,/x)-plane, this corresponds to the expected diagram Fig. [I] (left). A completely 
different scenario arises if instead the first-order region shrinks as /I is turned on, Fig. |3| (right). 
Now the physical point remains in the crossover region for any ju. Hence, existence and location 
of the critical point depend on two parameters of this phase diagram: the distance of the physical 
point from the critical line at /x = 0, and the sign and strength of the critical surface's bending. 

2.1 How to "measure" the properties of a phase transition 

A lattice calculation of the properties of a phase transition consists of two steps. First, one 
needs to find the phase boundary, i.e. the critical coupling /3 C separating the confined/chirally bro- 
ken regime from the deconfined/chirally symmetric regime. To this end one computes the expecta- 
tion value of the "order parameters" O for those symmetries, like the chiral condensate (yy) or the 
expectation value of the Polyakov loop, as well as its associated fluctuations, or generalized sus- 
ceptibilities, X = V((0 2 ) — (O) 2 ). For fixed quark masses, the transition is then announced by the 
rapid change of the order parameter as a function of lattice gauge coupling j8, and the peak of the 
susceptibility serves to locate it more precisely, cf. Fig. || (left). Using the two-loop beta-function 
(or some numerical fit to a non-perturbative beta-function), j8 c can be converted to a critical tem- 
perature T c . This first step of locating the phase boundary is comparatively easy. 

However, on a finite volume, this critical coupling, viz. temperature, is only pseudo-critical, 
no true non-analytic phase transition can occur on finite volumes. The order of the transition thus 
has to be determined by a finite size scaling analysis. With larger and larger volumes, the peak 
in the susceptibility will begin to diverge if there is a true phase transition, and the rate of its 
divergence, associated with some critical exponent, determines whether it is a first or second order 
phase transition. In the case of a crossover the peak stays finite and analytic. Another, for our 
purposes more practical, observable is the Binder cumulant, to be evaluated at j3 = j8 c , 



B A (m,n) = -^^ , 8 ¥¥ = ¥¥ -( ¥¥ ). (2.1) 



In the infinite volume limit, B4 — > 1 or 3 for a first order transition or crossover, respectively, 
whereas it assumes a particular value characteristic of the universality class at a critical point. 
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For 3d Ising universality, one has B4 — > 1.604. Hence B4 is a non-analytic step function, which 
gets smoothed out to an analytic curve on finite volumes, its slope increasing with volume to 
gradually approach the step function, Fig. |] (right). Once the universality class of a critical point 
is established by finite size scaling, B4 allows to investigate the change of the order of a transition 
already on a single (large enough) volume. In [ 13], we have used 8 3 , 12 3 and 16 3 x 4 lattices with 
standard staggered fermions to estimate the critical couplings for which B4 = 1 .604, and check their 
consistency on different volumes. Here, I also report on new jj. = simulations on 12 3 , 18 3 x 6. 



3. The chiral critical line m c s {m u ^) 

For Nf = 3, and at zero density, the critical quark mass marking the boundary between the 
crossover and the first order region is large enough for simulations to be carried out, and it was 
possible to determine it with some accuracy JT5| , [l(], [17|, 13]. Moreover, finite size analyses estab- 



lished it to belong to the Z(2) universality class of the 3d Ising model []15[]. The critical mass can 
be read off from the red squares in Fig. || (left), where data for the Binder cumulant are shown as a 



function of varying quark mass, anf = 0.0263(3) or m c /T c = 0.105(1) [13]. 

Starting from this point, the boundary line has been mapped for the non-degenerate Nf ■ = 2 + 1 
theory, again with N t = 4 [|13|], where we have assigned a chemical potential to the two degenerate 
light quarks only. The methodology is the same as in the three flavour case: fix a strange quark 
mass am s and scan the Binder cumulant in the light quark mass am u j for the corresponding critical 
point. Repeating for several strange quark masses produces the critical line am c s (am u j) shown in 
Fig. I (right). 

In order to set a physical scale for the boundary line, spectrum calculations at T ~ have been 



performed [13] at the parameters indicated by the arrows in Fig. H (right). They show that m s at 



the upper arrow is approximately tuned to its physical value (^ ~ ^| P hys=0.65), while the pion 
is lighter than in physical QCD (^ = 0.148(2) < ^| p h ys = 0.18). This confirms that the physical 
point lies to the right of the critical line, i.e. in the crossover region. 
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Figure 5: Left: B4 as a function of quark mass and /I, for Nf = 3, N t =4. Right: Comparison of the Binder 
cumulant for = on jV, = 4,6 lattices. 



4. Strong cut-off effects on the critical line 

While the universality class of a critical point is determined by long range fluctuations and 
thus should be insensitive to cut-off effects on a coarse lattice, the critical quark mass is a quantity 
requiring renormalization and expected to be highly sensitive. Here we present as yet unpublished 
data of a simulation on the finer N t = 6 lattice for Nf = 3. Fig. || (right) shows a comparison 
between the N t = 4 and N t = 6 data for B<n{m), as a function of m/T. If the lattice spacings were 
in the scaling region, the critical values extracted from the intersection with the Ising value should 
be close to each other. In contrast, the critical bare quark mass on the N t = 6 lattice is smaller by 
almost a factor of five, nf /T pa 0.02. We have also computed the pion mass corresponding to this 
nf on N t = 4,6, finding nf n rj 1.68,0.96, respectively. Similar observations have been made when 



comparing unimproved and improved staggered fermion results on N t = 4 lattices [15]. Hence, the 
boundary line in the phase diagram appears to be still far from the continuum result. Note that this 
implies the gap between the physical point and the critical line to widen significantly towards the 
continuum, thus moving He to larger values for fixed curvature of the critical surface. 



5. Finite density phase diagram using imaginary ji 

For imaginary /I = the fermion determinant is real positive and simulations can be carried 
out in the same way as at \i = 0. The problem then is to get back to the physical situation or real 
jj.. This is easily possible if an observable, measured as a function of imaginary chemical potential, 
can be fitted by a truncated Taylor expansion in fli/T, 

Note that, because of an exact reflection symmetry of QCD under /J. — ► — pi, only even terms appear. 
Clearly, /J. / (tiT) has to be small enough for this series to converge well, but there is a limited radius 
of convergence: when going in the imaginary direction, a Z(3) transition to the neighbouring, and 
unphysical, center-sector of the gauge group is crossed at fXj/T = 7i/3, thus limiting the useful 
range for this approach to \fi\/T<l. 
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5.1 The critical temperature at finite density 

The peak position in the susceptibilities calculated for several imaginary /I can be fitted by a 
leading 0{pL 2 ) and next-to-leading 0(jU 4 ) Taylor series, the data are well described by the leading 
term, fi c (am,an) = 5.1369(3) + 1.94(3) (am - am c Q ) + 0.781(7)(a/i) 2 + .. . Converting to contin- 
uum units by means of the two-loop beta-function yields the result 

Tc{ ^ m) -l+2.111(17)" c(0) -0.667(6)f I V + - (5.2) 



T c (n = 0,m c (0)) v ' nT c w \nT c (0,m) 

Note that one obtains coefficients of order one when expanding in the "natural" parameter (/x / kT), 
as one might expect on the grounds that the Matsubara modes of finite temperature field theory, 
coming in multiples of 7lT, are setting the scale of the problem. The coefficients depend on Nf, a 



comparison of the leading /^-coefficient for various cases can be found in [JlOj]. The curvature gets 



stronger with increasing Nf, which is consistent with expectations based on the large 7V C expansion. 



Subleading coefficients are also beginning to emerge as we discuss in section p.3[ Note that the 
continuum conversions using the two-loop beta-function are certainly not reliable for these coarse 
lattices, while fits to non-perturbative beta-functions tend to increase the curvature. Finally, detailed 
comparisons of /3 c (/l) calculated via reweighting techniques or analytic continuation have been 
made for Nf = 2,4 using equal lattice actions, and complete quantitative agreement is found [§]. 

5.2 The critical point as a function of quark masses 

Let us now move to the second step in the program and discuss the order of the phase transition 



as a function of quark masses and chemical potential. The strategy followed in [ |13| ] is to compute 
the change of the critical line as a function of imaginary chemical potential and fit it by a truncated 
Taylor series, just as was done for the critical coupling. Again the Binder cumulant was used to 
determine the order of the transition. The data for Nf = 3 and various /^ are shown in Fig. |5| (left). 
The chemical potential is found to have almost no influence on B4 as a function of quark mass. 
The change of nf with jx 2 , i.e. the leading coefficient in a Taylor expansion, is determined from 
d{anf) j d{a\x) 2 = -(dB^/ d{a}l) 2 ) / (dB^/ d{am)) and we obtain 

am c (a\i) =0.0270(5) -0.0024(160) (a\l) 2 . (5.3) 

Care must be taken for the conversion to physical units. The crucial point is that T = 1/ (aN t ), 
and for fixed N t the lattice spacing is adjusted with changing /i to tune the temperature so as to stay 
on the critical line, i.e. a(T(fi)) changes. Expressing the change of the critical quark mass with 
chemical potential in lattice and continuum units as 

anf(ii) 1 c, . m c (a) . / A 4 \ 2 ,r as 

am c (0) am c (0) m c (0) \%T J 

then c\ and c\ are related by 

r = %1 C J + ( 1 dT c {nf{n\n) \ 
1 Niamey \T c (m,il) d(^i/nT) 2 J ^ Q ^ } 

Since c\ is consistent with zero, the second term dominates, leading to an overall negative coeffi- 
cient c\ = —0.7(4) Jl3[]. This is evidence that, in the Nf = 3 theory on an N t = 4 lattice, the region 
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Figure 6: Finite difference calculation of the dfi c /d(ajj,) 2 (left) and c\ ~ (dB4./d(ajj,) 2 (right). The error 
bands correspond to the leading and next-to-leading order fits from Jl3[], 



of first-order transitions shrinks as a baryon chemical potential is turned on, and the unexpected 
scenario of Fig. |3| (right) is realized. Interestingly, similar qualitative conclusions are obtained 



from simulations of the three flavour theory with an isospin chemical potential [18], as well as 



simulations at imaginary /I employing Wilson fermions [19]. 



5.3 Check of systematics 

Clearly, this is a rather consequential finding, and one must ask how reliable it is. Firstly, 
our result for c\ is only two standard deviations away from one with opposite sign. Second, one 
might worry about systematic errors when fitting only the first term of the Taylor series. Higher 
terms with alternating signs could fake a leading coefficient c\ f» in a truncated series. Such terms 
would no longer cancel after continuation to real }i, leading to a different picture. In order to clarify 
this question we are currently performing a direct calculation of the derivatives d$ c jd(a\i) 2 and 
dB4/d(a/j,) 2 , thus avoiding a fit of the whole function by a truncated Taylor series. We compute 
the derivatives as finite differences, 

do _ Um oM-o(Q) ! (5 _ 6) 



d(afl) 2 (afi) 2 ^o (an) 2 

which offers several means to improve the signal. 0(a}i) is evaluated by reweighting, which has 
no overlap problem for \afx\ < 0.1, thus eliminating uncorrected fluctuations in the terms of the 
difference. The reweighting is done in the imaginary pL -direction, since then there is no sign prob- 
lem and the reweighting factor is less noisy. This in turn permits to compute the reweighting factor 
by noisy estimators rather than exactly, which is much cheaper. 

The results of these computations are shown in Fig. ^. The derivatives are obtained by extrap- 
olation to jit = 0: dp c /d(a/J.) 2 = —0.746(1), c[ = —0.11(2). Also shown in the figures are the one 
sigma error bands from our previous 0(jJ. 2 ) and next-to-leading 0(pL A ) direct fits to B4, as tabulated 
in [13]. Note that the derivatives are completely consistent with the \i A fits, but not with the leading 
order ones, i.e. we are indeed sensitive to more than just the leading coefficient. In the derivative 
calculation this is indicated by the slight slope. Note, that the errors of the derivative computation 
are dramatically smaller, so that now c\ is non-zero and negative beyond reasonable doubt. As 
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Figure 7: Left: Effect of keeping the quark mass fixed in lattice units in jTT|]. Right: Comparison at finite pt 
between the Nf = 2 + 1 and the Nf = 2 theory considered in fll^ ] . 

discussed above, the conversion to continuum units adds a dominant negative contribution to this. 
We conclude that, on an = 4 lattice, the phase diagram Fig. || (right) is realized. 



6. Discussion 

Our results about the critical surface appear to be in qualitative contradiction with those of 
jnj], [jl2|], which both conclude for the existence of a critical point (}JLe,Te) at small chemical 
potential He/Te < 1. However, in considering the reasons for such disagreement, one can see that 
the different data sets are actually not inconsistent with each other, and the differing pictures can 
be explained by standard systematic effects. 



In [ ]12| ] the critical point was inferred from an estimate of the radius of convergence of the 
Taylor expansion of the free energy. Regardless of the systematics when only 4 Taylor coefficients 
are available, the estimate is made for Nf = 2. The (ju, T) phase diagram of this theory might well 
be qualitatively different from that of Nf ■ = 2 + 1 QCD, as illustrated in Fig. [7] (right). Its critical 
endpoint point, obtained by intersecting a critical surface when going up vertically from the Nf = 2 
quark mass values, is clearly a long way from the critical endpoint of physical QCD, and nothing 
follows quantitatively from the value of one for the other. 

In [|TI|] the double reweighting approach was followed. By construction, this reweighting 
is performed at a quark mass fixed in lattice units: am uc j = N t ^- = const. Since the critical 
temperature T c decreases with jU, so does the quark mass. This decrease of the quark mass pushes 
the transition towards first order, which might be the reason why a critical point is found at small 
}X. This effect is illustrated in the sketch Fig. [7] (left), where the bent trajectory intersects the critical 
surface, while the vertical line of constant physics does not. Put another way, [ |Tl| ] measures the 
analogue of c\ instead of ci in Eq. (5_4). From their Fig. 1 [11] we see that the distance of the 



theory from criticality stays constant for small ju, consistent with a small or zero coefficient c\ « 0. 
Taking the variation a(r(ju)) into account could then make a dominant contribution, which might 
possibly change the results qualitatively. 

A robust finding is the high quark mass sensitivity of the critical point: irrespective of the 



sign, if C[ ~ 0(1) in Eq. (5.4), then m c (/i) is a slowly varying function of jj,, just as the pressure, 



screening masses or T c . Hence, /^(m) is rapidly varying. A change of quark masses by a few 
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percent will then imply a change of He by 0(100%), which makes an accurate determination of He 
a formidably difficult task. Our combined results of a negative curvature of the critical surface as 
well as an increasing gap between the /I = critical line and the physical point make the existence 
of a critical point at small He very unlikely. 

A last note of caution concerns the fact that most investigations have used unimproved stag- 
gered quarks on coarse N t =4 lattices only. This might be worrisome given the exceedingly light 
quarks involved, as one should take the continuum limit before the chiral limit for this discretization 



21]. Investigations on finer lattices are indispensable in order to settle these issues. Finally, 
a more complicated phase structure with additional critical surfaces in Fig. [| might allow for a 
critical point not continuously connected to the /I = transition. 

Acknowledgement: Unpublished work reported here is done in collaboration with P. de Forcrand. 
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